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Abstract 

We extend our Wilsonian renormalization group (RG) analysis on the pionless nuclear effective 
theory (NEFT) in the two-nucleon sector in two ways; on the one hand, (1) we enlarge the space 
of operators up to including those of 0{p'^) in the S waves, and, on the other hand, (2) we consider 
the RG flows in higher partial waves {P and D waves). In the larger space calculations, we find, 
in addition to nontrivial fixed points, two "fixed lines" and a "fixed surface" which are related to 
marginal operators. In the higher partial wave calculations, we find similar phase structures to 
that of the S waves, but there are two relevant directions in the P waves at the nontrivial fixed 
points and three in the D waves. We explain the physical meaning of the P-wave phase structure 
by explicitly calculating the low-energy scattering amplitude. We also discuss the relation between 
the Legendre flow equation which we employ and the RG equation by Birse, McGovern, and 
Richardson, and possible implementation of Power Divergence Subtraction (PDS) in higher partial 
waves. 
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I. INTRODUCTION 



Nuclear effective field theory (NEFT) is a low-energy effective field theory of 

nucleons based on the general principles of quantum field theory and the symmetries of the 
underlying theory of hadrons, QCD. (See Ref. , 0] for reviews.) At very low energies, 
where even the pions are regarded as "heavy," interactions of nonrelativistic nucleons are 
simulated by (infinitely many) contact operators. Such a theory is called pionless NEFT 
with the physical cutoff scale being around the pion mass, while the pionful NEFT is needed 
at higher energies, where the effects of the exchange of pions are explicitly taken into account. 
Because NEFT contains infinitely many operators, one needs an organizing principle, called 
power counting, to systematically calculate physical quantities to a certain order. 

It is interesting to note that the actual two-nucleon system in the S waves is fine-tuned. 
The scattering lengths are unnaturally large compared to the scale characteristic to the two- 
nucleon interaction. From the RG point of view, this unnaturalness may be rephrased as that 
the system is very close to the nontrivial fixed point (or better, to the critical surface) [3, Q]- 
The two-nucleon system is thus nonperturbative due not only to the strong coupling, but 
also to the closeness to the critical surface. 

The unnaturalness of NEFT makes the power counting issue complicated. The so-called 

h 

Naive Dimensional Analysis (NDA) |9j works for perturbative systems, but does not account 
for fine-tunirig. One needs a power coun ting which encodes the fine-tuning. There are a lot 
of papers 0, Isl, llO, 11, 12, 13, 14, 15, 16, 17\ devoted to the power counting issues in NEFT, 



and it is still an important subject of discussions. 



In a previous paper 



SIS 
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two of the present authors performed a Wilsonian RG analy- 



22| of pionless NEFT to determine the power counting of the operators in 



the ^5*0 and ^5*1-^-01 channels on the basis of the scaling dimensions at the nontrivial fixed 



point. We employed the Legendre flow equation [23|, |2j, |25|, |26|, |27|] as a formulation of 
nonperturbative RG, and we reproduced known results obtained by Birse et al. 28|] up to 
0{p'^) in a completely field theoretical fashion. We also emphasized the phase structure and 
identified the inverse of the scattering length as the order parameter. In the strong coupling 
phase, there is a coupling which grows as the floating cutoff A is lowered, and there is a 
bound state, while in the weak coupling phase, all the flows run into the trivial fixed point 
and there is no bound state. The determination of the power counting on the basis of the 
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scaling dimensions is a unifying principle of EFT, useful especially for the cases in which 
nonperturbative dynamics is important. 

One of the good features of Wilsonian RG analysis of EFT is that it provides an overview 
of all the possible theories consistent with the symmetries, not restricted only to the one 
that describes the real world. In other words, it characterizes the physical system in a broad 
perspective in terms of the RG flows. 

Although the formulation of the Legendre flow equation does not contain any approx- 
imation, in order to solve it, one needs an approximation; the restriction of the space of 
operators. It is important to note that the restriction does not require an a priori power 
counting more than just the canonical dimensional counting. The point is to include a 
sufficient number of operators. The scaling dimensions are determined for the linear combi- 
nations of the operators included at each fixed point. 

It is important to note that when the space of operators is enlarged, the results (e.g., 
the scaling dimensions) are expected to always converge to the true values. The question 
is how fast they converge. In order to see it, it is necessary to actually enlarge the space 
of operators. If the enlargement does not alter the results of the calculations with smaller 
space very much, we may conclude that the results are close to the true values. It has 
been argued that the small dependence on the choice of cutoff functions also suggests fast 



convergence 



29|. 



Since our formulation is quite general, we can readily do a similar analysis to other sys- 
tems. For example, one can investigate higher partial waves, where higher derivative terms 
play an important role. Note however that, in general, higher partial waves are physically 
not so significant at very low energies and that the physical two-nucleon system does not 
seem to be fine-tuned in those higher partial waves. Nevertheless, it is interesting to perform 
the RG analysis for higher partial waves, because it provides a better characterization of 
the physical two-nucleon systems than just that of the angular momentum. (How many 
nontrivial fixed points there are? How far is the physical system from them?) This kind of 
information may be useful in the study of other systems described by a similar EFT. 

One may think of the so-called "halo nuclei" as such an example. A typical halo nucleus 
consists of a core (e.g., alpha^particle) and a few "halo" nucleons, and an EFT is proposed 



to describe such a system 



30|. A particularly interesting example is the nucleon-alpha 



scattering in the P-wave, for which the P3/2 channel displays a resonance just above the 
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threshold 



30 



311 ] ■ Such a resonance is clearly related to a RG critical surface. The RG 



structure for nucleon-alpha scattering is very similar to that of the pionless NN scattering 



which we consider in this paper. We provide an account on the power counting of Refs. 30|, 



3l[ | on the basis of Wilsonian RG analysis 

In this paper, we continue our previous study and perform a Wilsonian RG analysis 
of pionless NEFT in two-nucleon sector. After recapitulating the previous paper, we first 
enlarge the space of operators to include those of 0{p'^) in the ^5*0 channel. We find two "fixed 
lines" and a "fixed surface" besides trivial and nontrivial fixed points. The nontrivial fixed 
point that we think the most relevant to the real world persists, and the scaling dimensions 
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of the eigen-operators of the linearized RGE do not change, consistent with Ref. 

Second, we consider the P waves and obtain the phase structure by solving the RG 
equations. The phase structure is similar to that in the 5* waves, but there are two relevant 
operators in the P waves at the nontrivial fixed points. We also discuss briefly the D waves 
and find that there are three relevant ones. 

The structure of the paper is the following: In Sec. [Tll we recapitulate the main points 
of the previous paper 18|] in order to introduce the notations and the main concepts in our 
a„alys,s. Smce the present paper closely followsthe d,scu.s,ons given m Ref. Q. we expect 
that the readers intimately consult with Ref. [18|] . In Sec. IIIII we discuss the enlargement of 
the space of operators. The calculations for the higher partial waves are given in Sec. IIV[ 
We also discuss possible extensions of Power Divergence Subtraction (PDS) renormalization 
for higher partial waves. The summary and discussions are given in Sec. |Vl In Appendix 
we discuss the relation between the Legendre flow equation in the sharp cutoff limit and 



the RG equation used by Birse et al. 
studied in Appendix [Bl 



28 



. The cutoff function dependence of the results is 



II. NONPERTURBATIVE RG ANALYSIS FOR THE TWO-NUCLEON SYSTEM 
AT VERY LOW ENERGIES 

A. Power counting and the scaling dimensions 

The most basic idea behind the power counting is the order of magnitude estimate based 
on dimensional analysis. For the two-nucleon system at very low energies described by the 
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pionless NEFT, there is a physical cutoff scale Aq ~ 0{m.^), above which the effective theory 
is not applicable. Since it is the only scale (except for the nucleon mass, M, see the footnote 
3 of Ref. 1^ and the section 1.2.1 of Ref. 32]) which appears in the pionless NEFT, the 
dimensional analysis for the pionless NEFT is based on this scale. 

Quantum fluctuations may change the classical dimensional analysis. The quantum 
counter part of the dimension is the scaling dimension, which can be obtained by the RG 
analysis. It is therefore natural to consider the power counting based on the scaling di- 
mensions. Wilsonian, or nonperturbative, RG is a suitable tool to handle the quantum 
fluctuations. 

It is known that there is another scale in the two-nucleon system, the scattering length 
in an S wave, the inverse of which is much smaller than Aq. It can be understood that such 
a scale is not fundamental, but is a result of the fine-tuning of the parameters of the EFT. 
In the RG language, the fine-tuning is closely related to the existence of a nontrivial fixed 
point and a critical surface of the RGE. It has been shown that on the critical surface the 
scattering length is infinite. In order to get a large scattering length, the coupling constants 
must be fine-tuned to be near the critical surface. 

Around the nontrivial fixed points the scaling dimensions are drastically different from 
their classical values, the canonical dimensions. The coupling which corresponds to the 
scattering length classically has negative dimension, indicating that it gets less important at 
lower energies. But near the nontrivial fixed point (or the critical surface), the (quantum) 
scaling dimension is positive, so the corresponding interaction becomes more important at 
lower energies. On the other hand, near the trivial fixed point, where all the coupling 
constants are small, the scaling dimensions are the same as the classical values. There the 
quantum fluctuations do not alter the importance of the interaction. 

A relevant operator, which has a positive scaling dimension, controls the deviation of the 
coupling constants from a critical surface. In order to get a scale much smaller than Aq, 
the relevant operator must be fine-tuned. Near a nontrivial fixed point, the number of the 
parameters to be fine-tuned is that of the relevant operators. 

In determining the power eounting, one usually does two things: (1) choose a set of 
operators one is going to work with, and (2) determine the importance of the operators. In 
a conventional power counting scheme, one first determines the importance of the operators 
(2) by assuming some scaling property, and it then leads to the set of operators (1) to 
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a given order. On the other hand, in the Wilsonian RG approach we only choose a set of 
operators (1). The importance of the operators (2) is determined as an output by the scahng 
dimensions of them. 

In most cases, it is unhkely to miss important operators in both approaches, as far as one 
includes all the operators of low canonical dimensions to a certain order. What matters is 
the determination of the importance of the operators (the scaling property). 

Up to the point where the fixed points and the scaling dimensions are obtained in the 
Wilsonian RG approach, however, no contact has been made with where the real world is in 
the RG flow. The Wilsonian RG approach reveals all the possible theories described by the 
same action without the knowledge of the actual systems. It is the strength as well as the 
weakness of the approach. Namely, one can determine possible scaling properties consistent 
with the action (the symmetry, spacetime dimensions, and the degrees of freedom) while 
one still needs additional information about the system (e.g., a large scattering length) one 
is trying to describe, in order to determine where it is in the RG flow. 

Convergence of the approximation has different meanings in the conventional and the 
Wilsonian RG approaches: in the conventional approach, the validity of the power counting 
is examined (in most cases, numerically) whether the expansion converges. With a wrong 
assumption of scaling property, the expansion fails to converge, i.e., even though one includes 
more operators, the results (the phase shift of a specific channel, for example) do not improve. 
See Ref. [isl for the important example. On the other hand, in the Wilsonian RG approach, 
the results (the scaling dimensions at a specific fixed point, for example) improve as the 
space of the operators is enlarged. The question is how fast the results converge. 



B. Legendre flow equation 



Legendre flow equation [23 



24 



25l . 1261 . l27j is one of the implementations of Wilsonian 
RGE. It is formulated as a RGE for the infrared (IR) cutoff effective action rA[$] called 
effective averaged action, in which the quantum fluctuations above the cutoff have been 
integrated. The averaged action is the generator of the one-particle irreducible (IPI) vertex 
functions containing only the fluctuations p > A. 
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The Legendre flow equation is given by 



dA 



Ilk 



(2.1) 



where r(2) stands for the second derivative of the averaged action Fa with respect to the 
fields, and is the inverse of the full propagator dressed by the vertices with classical field 
insertions, containing only the fluctuations p > A, and the Tr denotes the integration over 
momentum and also the trace in the internal space. The function i?A effectively cuts off the 
IR part of the fluctuations. Our choice is, as in the previous paper, as follows; 



Ra{p' 



PL 

2M 



exp 



El 

A2 



(2.2) 



where p is the three-momentum. In the n — oo limit, it becomes a sharp cutoff. Another 



choice is adopted in Ref. [33|. We derive the RGE for an arbitrary value of n, but the results 
look so complicated that we present only the results in the n —>■ oo limit. The n- dependence 
is studied in Appendix [Bl 

In the application to the two-nucleon system, we have drastic simplification of the RGE. 
(i) Because of the nonrelativistic feature of the system, we do not include anti-particles. 
Thus there are no self-energy corrections, nor tadpole contributions, (ii) In the two-nucleon 
sector only the four-nucleon (4N) operators contribute. From these, we end up with the 
one-loop diagrams involving two vertices with the tree-level propagators contributing to the 
Legendre flow equation. The multi-loop effects are encoded in the cutoff dependence of the 
(infinitely many) coupling constants. 

In Appendix [Aj we discuss the relation between the Legendre flow equation in the sharp 



cutoff limit and the RG equation used by Birse et al. for the potential V, 



dV 
dX 



M 



V{k\A,p;A)- 



A2 



;V{A,k,p;A). 



(See Ref. 



27r2 ' ' 'A2-p2 
for details.) It is shown that they are essentially equivalent. 



(2.3) 



As we emphasized in Ref. [18[, there is no obvious way of imposing a Galilean invariant 
cutoff at the averaged action level 331], so that the inclusion of the IR cutoff function in the 
averaged action should be understood as a symbolic one. Note that the IR cutoff function 
constrains the momenta of individual particles, but they are not invariants under Galilean 
transformations. Fortunately, however, the correct way of implementing a cutoff is clear in 
the two-nucleon system; to impose the cutoff on the relative three momentum. See Sec. 4.1 



in Ref. [18j for the explicit manipulation. Note that it is not an approximation, but the only 
way that we know to implement the cutoff in a Galilean invariant way. 



C. Pionless NEFT up to ©(p^) in the ^Sq channel 



Although the Legendre flow equation is exact, one needs an approximation to solve it. 
We consider a simple truncation of the space of operators. We retain only the operators with 
derivatives up to a certain order. We simply count the number of spatial derivatives (V ~ p) 
and a time derivative is counted as two spatial derivatives {dt ~ p^). The approximation is 
based on our hope that, even though some operators get large anomalous dimensions, their 
"ordering" of importance would not change very much; the lower the canonical dimension 
is, the lower the scaling dimension would be. We also expect that there are mixings among 
operators. In Ref. 18|, we consider the following ansatz for the averaged action up to 0{p'^), 



+ C2 
+ 2B 



t 



jYTpe5o)jYj (^A^^pj'^o) V^iv) + h.c. 



(2.4) 



for the ^5*0 channel, where Pa ° is the projection operator to the channel, 

pi'So) ^ 2 a 

V8 ' 



(2.5) 



and V ^ = + — 2 V ■ V. As we emphasized in a previous paper [3J], it is important 
to include the so-called "redundant operators," which can be eliminated from the action by 
using the equations of motion. 

By introducing the following dimensionless coupling constants. 



MA^ MA\^ A3 



27r2 



27r2 



27r2 



(2.6) 
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the RGEs can be expressed as 



dx 

Iti 

dy 
dt 

dz 
'dt 



—X — 



-3y 



-3;^ + 



+ 2xy + y"^ + 2xz + 2yz + z^ 



+ 2xy + ^yz- ^z^ 



^x^ + xy + ^y^ -xz-yz- 



(2.7) 
(2.8) 
(2.9) 



where t = In (Aq/A). We found a nontrivial fixed point, 

1 1 



(x*, y\ z*) 



' 2'2 



(2.10) 



as well as the trivial one (0, 0, 0). The eigenvalues and the corresponding eigenvectors of the 
linearized RGE at the nontrivial fixed points are found to be 



z/i = +l: Wi = (1,1,-1), 
z/2 = -l: = (0, -1, 1) , 

z/3 = -2: = (2, -1,-2). 



(2.11) 
(2.12) 
(2.13) 



(Note that the signs and the normalizations of the eigenvectors are arbitrary.) The eigen- 
vector associated with the positive eigenvalue ui corresponds to the scattering length. 
In Ref. 3J] we performed a similar RG analysis based on the cutoff independence of 



the amplitude and found the same scaling dimensions with slightly different eigenvectors. 
Since eigenvectors are not universal quantities and the approximations are different, it is 
therefore not surprising that the eigenvectors are different. The important thing is that the 
eigenvalues, which are considered as the universal ones, actually agree. The fact that the 
first two eigenvectors agree refiects that these two approximations are similar. 

These results are obtained with a severe restriction of the space of operators. One should 
examine the validity of the approximation by actually enlarging the space of operators. It 
is expected that some properties, i.e., the scaling dimensions at the nontrivial fixed points, 
are universal and would not change very much under the enlargement if the approximation 
is good, while the others such as the directions of eigenvectors may change. It is well known 
that the restriction of the space of operators causes various artefacts 35, [s^l- 

In the next Section, we investigate the effects of the enlargement of the space of operators 
and compare the results with the previous ones presented in this section. 



III. PIONLESS NEFT UP TO ©(p^) IN THE CHANNEL 



A. Independent operators 

In this section, we enlarge the space of operators in the ^5*0 channel up to including those 
of 0{p'^). (The results for the ^Si-^Di channel are essentially the same as those for the ^5*0, 
but a bit more complicated.) We consider the following ansatz for the averaged action, 



rip*) _ r(p') 

J- A — ^ A 



+ / d^X 



-a 



41 



C42 (iV^Pi'^^^V^iv)^ (^N^Pp^^V^N 



2Bi 



2M 



N )■ +h.c. 



2M 



N 



N } +h.c. 



(3.1) 



We emphasize that the ansatz given above contains all the independent operators consis- 
tent with Galilean invariance, parity, spin and isospin invariance to the given order. There 
are actually two other operators which satisfies these requirements. 



N 



h.c. 



2M 



h.c. 



(3.2) 
(3.3) 



but they are written as linear combinations of the operators contained in the averaged action 
up to total derivatives, and thus we dropped them. The relations may be expressed most 
clearly in momentum space. Using momentum conservation, (pi + = (ps + P^Yi we 
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obtain 



ri2 {Si + S2) + (5*3 + S4) =ri2 {Si + S^j + r^^ {Si + S2) 



(^?2+ ^34) +7777^12^34, (3.4) 



AM ^ 2M 

4 



S1S2 + 5'3S'4 — (5*1 + 5*2) (5*3 + 5*4) — - Si 



2 

i=l 



where we have introduced the notations, 

r,, = (p,-p,)% S,^p^-^. (3.6) 

It is also important to note that we have included the interaction {Bi) which depends not 
only on the total energy of the two nucleons, but on the individual energies. The potential 



corresponding to this interaction is not considered in Ref. 5 



We do not include the relativistic correction terms such as A^^ (V^/8M^) in the averaged 
action because it breaks Galilean invariance. It means that we are considering the system 
at very low energies so that the relativistic corrections can be neglected. The effect of this 
term is estimated to be smaller than any of the terms we have included in Eq. (13.11) due to 
the additional (9((A/M)^) suppression. 



B. RG equations 

which 



The Legendre flow equation may be obtained in a similar manner as in Ref. 
generahzes the set of RGEs ([22D, (123]), and fl^ . 

By introducing the following dimensionless coupling constants, 

MA^ A^ 
«^ = ^16C4.(^ = 1,2), ^^ = ^^^^(^ = 1'2), ^3 = ^4^3, (3.7) 

together with x, y, and z defined in (12. 6p . we have a set of the RGEs for these variables in 
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the n — s> oo limit, 



dx 
'dt 



dzi 



dZ3 

dt 



- X - x{x + 2y + 2z + 2ui - 2zi) - y{y + 2z + 2ui - 2zi) 

- z{z + 2ui — 2zi) — ui{ui — 2zi) — zl, (3.8) 

dy fl \ /3 

— = - 3y - X \-x + 2y + ui + U2 + Zl + Z3 j - y \-y + z + 2ui + U2 + Z3 

+ Z - U2 - 2Zi - 23^ - Ui +U2 + Zi + 23^ + U2Z1 

+ Zl + 23^ , (3.9) 

dz /I \ /I 

— = - Sz + X \-x + y - z + ui + Zl + Z2 - Z3 j + y i-y - z + ui + Zl + Z2 - Z3 

-.(|. + ..-3.,-.. + .3)+«.(i.. + .. + ..-.. 

- Zl(^Zi + Z2- Zs^ (3.10) 

dui 
"dt ~ 

^ = - 5m2 - x{x + % + 2ui + 2U2) - y{^y + 4mi + 4^2 + 2zi + 223) + 2(221 + 223) 
dt 



2 

— 5ui — ui{x + y + z + ui — 2i), (3-11) 



- Mi(mi + 2^2) - ^2(^2 + 22i + 223) - 2i(32i + 4^3) - zj, (3.12) 

- 52i — 2i(a; + y + 2 + ui — 21), (3.13) 

- 522 + x{x + 2y -2z + 2ui - 223) + y{y ~ 2z + 2ui - 223) 
+ z{z - 2mi - 421 - 422 + 223) + ui{ui - 223) + 2;i(32i + 42;2 - 223) 
+ 2:2(^2 -2^3) + ^3'. (3-14) 

- 523 + x{x + 3y - z + 2ui + U2 - 23) + y{2y - 2z + 3ui + U2 + 21 + 22 - 223) 

- 2;(mi + U2 + 3zi + 22 + 22:3) + ni(ni + U2 - Z3) + ^2(21 + 2:2 - 2:3) 

+ 2i(32i + 22:2 + 23) + 2223 - 2:3. (3.15) 

By setting ui = 1x2 = 21 = 22 = 23 = by hand in the first three of these RGEs and 
ignoring the rest, they reduce to the ones, (12.71) . (12.81) . and (12. 9p . as they should. This is 
exactly what we did in the O(p^) calculations. 

To analyze the RG flows, it is useful to determine the fixed points, at which the coupling 
constants do not run. We use Mathematica to solve the fixed point equations and found the 
following solutions. 
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23125-1000Mi-1700y+36y^ on_L,, 19375- IOOOmi +300y- 36y^ -1875+1000};-36y^ 



A (0,0,0,0,0,0,0,0), 
R ('—1 — i i n — ^ ^ i\ 

^^ 1, 2' 2' ^' 3' ' 3' 3/' 

C 123 21 Q 92889 n 13689 33831 \ 

^ V 4' 44 '22'^' 15488' ' 15488' 15488/' 

D (-25,^,-f,«i,20-2«i,-l^,MHi,H)) 

E (O, 0, 0, 0, ^2, 0, X^, -5X^+Xl+,u,+2X^u,+u,- 
F (-2'^ y -y y, 23125-1000^.1-1700^+36^^ -20 + ?yi 

where we have introduced 

X±=5±2V5y/^-U2. (3.16) 

The first fixed point [A] is the trivial fixed point. The second one [B] corresponds to the 
fixed point we considered in the previous paper (the fixed point Eq. fl2.10p ) because their 
first three coordinates coincide with each other. Note that, since the restriction of the space 
is nothing but the projection to the lower dimensional subspace, the fixed point Eq. fl2.10p 
has the components which we naturally expect the corresponding fixed point to have if the 



approximation is good. As we explained in Ref. [18[, this fixed point is the most relevant 
to the real two-nucleon systems, and we will discuss it later in detail. The third one [C] 
is a strange fixed point. We calculate the scaling dimensions at this point, and found that 
they are: —8, —13/2, —13/2, 5, 3, 1, and (—5 ±2-\/4T)/2. Because of the appearance of the 
complex scaling dimensions, we suspect that this is an artefact of the truncation. We will 
give the argument shortly. 

The fourth [D] and the fifth [E] are actually "fixed lines." The [D] is a fixed point for an 
arbitrary Ui, and the [E] for an arbitrary U2 < 0. The [D] is a straight line, while the [E] is a 
parabolic curve with the peak at (0, 0, 0, 0, 0, 0, 5, 0). The sixth [F] is an even stranger "fixed 
surface." Interestingly the scaling dimensions are the same on the whole fixed lines (surface). 



The scaling dimensions at [D] are: 5,0, (ll±iV239)/6, -7.86665±3.58098i, and -1.35557± 
0.3096532. Since this also has complex scaling dimensions, we suspect that it is an artefact 
too. See below. The [E] has the scaling dimensions, — 5, — 5, — 5, — 3, — 1, 0, 2, 5. Note that 
all the scaling dimensions are integers, and the existence of a marginal direction. Actually, 



a fixed point on it has been noticed by Birse 



561 ]. The [F] has the scaling dimensions, 
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—5, (—11 ± v4T)/2, 0, 0, 2, 5, 7. Note that there are two marginal directions, which are 
tangential to the surface. 

We consider the fixed point [C] and and the fixed line [D], which have complex scaling 
dimensions, artefacts of the truncation. In the following, we argue for this statement by 
giving a list of observations. We admit that each of them is not completely convincing, but 
they all suggest that fixed objects with complex scaling dimensions are artefacts. Of course, 
there is no definite criterion by which we can decide whether a fixed object is an artefact. 

The observations are the foUowings: (1) It is well known that the truncation generates 



complex scaling dimensions 



361], and it seems a general feature of truncation. So the appear- 



ance of complex scaling dimensions can be seen as a signal of the artefact of the truncation 
and does not seem to have real physical relevance, at least in most cases. (2) In the litera- 
ture, it is extremely rare for complex scaling dimensions that do not seem to be artefacts to 
appear in physical systems. To our best knowledge, only hierarchical Ising models 



and 



gravitational collapse 
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40| are such systems. The former are frustrated at every length 



scale, while the latter is known to have critical limit cycle (discrete self-similarity). We see 
that these examples are very different from the simple system of self-interacting nonrela- 
tivistic fermions. (3) In a previous paper [l^ we found a fixed point with complex scaling 
dimensions, but when enlarging the space of operators in the present paper, we do not find 
the corresponding fixed point. It disappears. This is a concrete example of a nontrivial 
fixed point which appears with complex scaling dimensions in smaller space of operators, 
and is in fact an artefact of the truncation. It is important to note that the very existence 
of a nontrivial fixed point can be an artefact. (4) It seems unlikely to get complex scaling 
dimensions by the method employed by Birse et al. 28[ , which is proved in Appendix |A] to 
be essentially equivalent to ours. Note that, even though the RGE for the potential to be 
solved is shown to be the same for the both approaches, the way of solving it is different. 
In the approach by Birse et al. the scaling dimensions are the sums of the powers of the 
monomial perturbation around the (inverse of the) fixed point potential that is not expanded 
in powers of the energy. The scaling dimensions obtained in such a way should be integers. 
(5) Because of the nonrelativistic feature of the present theory, only power divergences are 
involved. It is therefore natural to expect integer scaling dimensions. From this point of 
view, the fractional scaling dimensions of the fixed surface [F] also seem to be artefacts. 
Again, we cannot completely exclude the logical possibility that the fixed objects persist 
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in further enlargement of space of operators, while their scaling dimensions become real. 

Even though their physical relevance is unclear at this moment, it is interesting to find 
such possibilities of several kinds of fixed objects appear in the system. These possibilities 



28 



for 



have never been revealed in other papers. Compare with the work by Birse et al 
example. They started with an ansatz for a fixed potential that depends only on the total 
energy, found the solution, and analyzed perturbations around it. In our analysis, on the 
other hand, we scanned the whole (though restricted) theory space and found all possible 
fixed objects. 

In any case, irrespective to whether these fixed objects are artefacts or not, the most 
relevant fixed point to the realistic two-nucleon systems is the fixed point [B]. It is easy to 



28| up to including 



Z/i = 


+1 : 




Z/2 = 


-1 : 


U2 


1^3 = 


-2 : 




1/4 = 


-3 : 




1^5 = 


-4 : 


tt5 


= 


-4 : 


Uq 


U7 = 


-4 : 


Uj 


Z/g = 


-5 : 


Us 



see that our solution is the same as the one obtained by Birse et al. 

At the nontrivial fixed point [B] , we find the following scaling dimensions and correspond- 
ing eigenvectors, 

11 11 11\ 

l,l,-l,0,y,0,-y,-yj, (3.17) 

(0,-1,1,0,-4,0,4,4), (3.18) 

^•^■4°-f-«--f-f)- 

(0,0,0,0,1,0,-1,-1), (3.20) 

1,-1,1,^,0,3,0,0^, (3.21) 

1,-1,1,-^6,0,6,0^, (3.22) 

:-2, -1,1, 5, -6, 0,0, 3), (3.23) 

1 1 1 22 4\ , , 

2, 0, -, 0, -,- . (3.24) 



2' 2' ' ' 6' ' 3 ' 3^ 

Note that (i) all the scaling dimensions are integers, (ii) The first three eigenvalues coincide 
with those obtained in the previous paper, which justifies the identification of the fixed 
point [B] in the present paper with the fixed point Eq. fl2.10p in the previous paper, (iii) 
The results are insensitive to the choice of the cutoff function. We give the analysis of the 
n-dependence in Appendix [Bl (iv) Only the eigenvector contains the zi component, 
i.e., depends on individual energies of two nucleons. This eigenvector is not considered in 
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Ref. [28 1 . 

In this section, the enlargement of the space of operators in the two-nucleon system in 
the ^Sq channel is considered. We found a new kind of extended fixed objects which have (a) 
marginal direction(s). We then concentrated on the fixed point which is the most relevant 
to the realistic two-nucleon system, and found that the universal properties we got in the 
previous paper are unchanged under the enlargement of the space of operators, confirming 
that the approximation is reliable. 



IV. HIGHER PARTIAL WAVES 

Operators which contribute in higher partial waves contain more derivatives, and thus 
are less important at low energies. In a few instances, however, higher partial waves contain 
some interesting information. A well known example is the P3/2 wave of the n-a system, 
where a narrow resonance state exists near the threshold, which is discussed in the NEFT 



context in Refs. 



30 



3l| . An important feature of this system is that there are two coupling 



constants to be fine-tuned. In the following, we will show that it comes out very naturally 
from our Wilsonian RG analysis. A similar analysis can be done for the D waves, and we 
show that, if there is a bound (or a resonance) state near the threshold, there are also three 
couplings to be fine-tuned. 

A. P waves 



In the NN system there are four channels in the P waves: ^Pi, ^Pq, ^Pi, and ^P2-^F2, 
where we consider the mixing with an F wave for the J = 2 channel. For the ^Pi channel, 
the interaction terms for the Fa may be written as 



-pint 
J- A 



C 



+ 2B, 



(4.1) 
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up to including O(p^), where is the projection operator to the ^Pi channel, defined as 



4V2 



Similarly, one can define the other projection operators. 



"\/3 

-^Qfci V k{io-20-i){ir2ra 

o 



,(3pi) 



ija 



8^2 



32 



1 



2 

V i(zcr2crj) + V j{ia2iJi) - -(5^^ V kii(J2crk) 



[lT2Ta 



V.V,V,-5^V.*,+ V/,.+ VA,)V 



They are normalized in the same way as in Ref. 



pol.avg 



13| 



{ia2(7i){iT2Ta) 



Tr 


'p{s) p(s')] 


=ilfel 






2 ' ' 



\2l 



(4.2) 

(4.3) 
(4.4) 
(4.5) 
(4.6) 

(4.7) 



where k is the relative three-momentum and / is the orbital angular momentum of the 
channel. For the P waves / = 1. The operators in another channel are obtained by replacing 



P/ with the corresponding projection operator. In the ^P2-'^P2 channel, there is an 



additional operator. 



I] a 



■I J a 



(4.8) 



which represents the mixing. Since the calculations are completely parallel, we will only 

n p \ 

demonstrate the results in the Pj channel in the following. 



We introduce the dimensionless coupling constants. 



2^2 -2 , y{l,l) - 27^2 ^4 , 

to write the following RG equations in the n ^ oo limit. 



^(1,1) 



2n 



2^4 



(4.9) 



dt 



dy{i,i) 
dt 

dt 



- 3x(i,i) - + 2x(ij)y(i,i) + 2x(i,i)Z(i,i) + y^i^i) + 2y(i,i)Z(i,i) + z^i^i)] 



5^(1,1) 



+ 2x(i,i)|/(i,i) + + Z/(i,i)2(i,i) - ^4,1) 



(4.10) 
(4.11) 



(4.12) 
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We obtain the RGEs of the same form for the ^Pq and ^Pi channels with the couphng 
constants being suitably defined. 

Compare them with the RGEs for the channel, ([22D, (ESD, and (^M>- With the 
appropriate replacement of the coupling constants, the coefficients of the quadratic terms 
are the same. Only the coefficients of the linear terms, which are nothing but the canon- 
ical dimensions of the corresponding operators, are different. This is because the channel 
dependence enters into the averaged action only through the projection operator and the 
canonical dimensions of the coupling constants. See the interaction part of Eq. (12. 4p and 
Eq. (14.11) . for example. 

One can easily obtain the fixed points, 

25 35 5^ 

'y'"6"'~2. 

We find that the scaling dimensions at the third fixed point are complex and may be disre- 
garded. 

It is useful to define the following variables. 



{x\y*,z*) = (0,0,0), 



9 _9 

'^'2' 2 



(4.13) 



u 



^(1,1)' 



2 (^(1.1) 



^(1,1)) 



w 



2 (^(1.1) 



^(1,1)) 



in terms of which the RG equations are written as 



du 

— = —6U 

dt 



[u + 2w)^ 



dv 
dw 



-5v — -(n + 4v){u 



2w) 



—— = —5w — w(u + 2w) 
dt 



(4.14) 

(4.15) 
(4.16) 
(4.17) 



The nontrivial fixed point is now at {u*, v*, w*) = (—3, |, O). Note that fiows starting in the 
w = plane never depart from it. 

In Fig. [U, we show the RG fiow in the w = plane. It is easy to see that u = —3 is 
a phase boundary; the fiows in the right of it (the weak coupling phase) run to the trivial 
fixed point, while those in the left of it (the strong coupling phase) go to infinity. 

At the nontrivial fixed point, we can easily obtain the following scaling dimensions and 
corresponding eigenvectors (in the u,v,w basis). 



^1 



3 : Ui 



-3 

V ° / 



^2 



1 : U2 



^3 



-2: W3 



12 
1 



(4.18) 
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RG flow for the P-wave 




u 



FIG. 1: The RG flow for the P wave in the w = plane. 



Note that there are two relevant directions, though they are all irrelevant at the trivial fixed 
point. For the first two eigenvectors, the scaling dimensions shift from their canonical values 
by six. 

We can investigate the ^P2-^-^2 channel in a similar way. After introducing dimensionless 
couplings, 

,5 



^(3,2) 



^(3,2) 



27r2 
AS 



2 ) 1/(3,2) 



2tx 



4MA5 (3p,) 
" 27r2 ^4 ' 
v^MA5^(3p,^) 



2 27r2 



4 ) 



(4.19) 



we have (in the n ^ oo limit) 

(^2^(3,2) _ 



dt 

dyi3,2) 
dt 

dZ(3^2) 

dt 

C?W(3.2) 

dt 



3X(3,2)- [a;j3,2)+ 2X(3,2)l/(3,2)+ 2X(3,2)^{3,2)+ 1/^3,2) + 2?/(3,2) ^(3,2) + 4,2) + 2m^?3,2)] , 

(4.20) 



- 52/(3,2) - 
-5^(3,2) + 



^a;(3,2) + 2X(3,2)2/{3,2) + ^y?3,2) + 2/(3,2) 2;{3,2) - ^4,2) + ^^3,2) 



, (4.21) 



-X 



(3,2) + 3^(3,2)1/(3,2) - a:(3,2)2(3,2) + -2/(3,2) " 2/(3,2)^(3,2) " -2^(3,2) + «^(3,2) 



- 5w;(3,2) - [a;(3,2)U'(3,2) + 2/(3,2)^^(3,2) + ^(3,2)^^(3,2)] 



(4.22) 
(4.23) 
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Note that, with the appropriate replacement of the couphng constants, the coefficients of the 
quadratic terms of the above RGEs are the same as those in the ^Si-^Di channel, and the 
only difference is the coefficients of the linear terms representing the canonical dimensions, 
just as we explained for the ^Pi and ^5*0 channels. The flows and the scaling dimensions are 
similar to those for other P waves, just as those for the ^Si-^Di channel to those of ^5*0. 



B. Amplitude for the P waves 



Similar results may be obtained by explicitly calculating the scattering amplitude as we 
did for the 5- waves in Refs. which is the extention of the method in Refs. [UJ, |4l|] 



to include the redundant operators. Consider the Lippmann-Schwinger equation with the 
"potential" in the center-of-mass frame (See Appendix R|) . 



-iV = -iP'^ ® P (-4pi ■ P2) 



C2 + 4C4(p?+p^)-2fi4(/ 



2M 



(4.24) 



where pi and p2 are momenta of the nucleons in the initial and final states respectively, P 
stands for the spin-isospin factor of the projection operator (for example, P = ^^0-2X2 for 
the ^Pi channel), and the ansatz for the amplitude, 

-zA= -iP^ ® P (-4pi ■ P2) [X(/) + F(/) {p\ + pI) + Z(/)p?p2] ^ (4,25) 

where X(p''), l^(p°), and Z{p^) are functions to be determined. The Lippmann-Schwinger 
equation reduces to 



- z^(p°,pi,p2) = -iV{p^,pi,p2) 



(2^ 



-iV{p^,k,p2] 



p^-\f+ie 



-iA{p\pi,k) 



where we have introduced 

A = [Xip') + Yip') {pi + pI) + Z{pViPI] 



V 



C2 + 4C,{pI+pI)-2bJp' 



,0 Pt+P2 

2M 



The solution is easily obtained as. 



A- = i 

D 



(C2 - 254/) + (4C4 + ^) /3 



Y 
Z 



D 



4Ca 



54 
M 



AC. 



M 



(4.26) 

(4.27) 
(4.28) 

(4.29) 
(4.30) 
(4.31) 
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with 



D=l- (C2 - 254/) h - 2 (^4C4 + h 



4 \ j2 



where 



M 



B, 
dk 



4C4 + §) /./3, 



/i 



^-Mp^ - ie. 



(4.32) 



(4.33) 



The RG equation may be obtained by requiring the amphtude to be independent of the 
cutoff A. By introducing 



2 t: 

X = l + -w, y = u- -w^, Z = 2{v + w) + -w^ 
5 



7 



we find the following RG equations, 



= (lOA" - + 2xy - y) , 

A— = ^ (3^' - SA'^Z + IQX^Z + 2A'3^Z) 
dA ^ ^ 

which give rise to the following fixed points, 

(A■^r,^*) = (1,0,0), (1,-3,9). 

In the original variables, they are 



(M^y^u;*) = (0,0,0), 



-3,^,0,, 



(4.34) 

(4.35) 
(4.36) 
(4.37) 

(4.38) 

(4.39) 



consistent with the previous analysis using the Legendre flow equation. Note that there are 
no additional nontrivial fixed points. 

At the nontrivial fixed point, we can obtain the scaling dimensions and the corresponding 
eigenvectors (in the m, f , w basis), 

1 1\ fo\ /l2^ 



Ui = 3 : Ui 



1^2 = 1 : U2 



V ° / 



1/3 = -2 : xt3 



-13 



(4.40) 



Compare the results with (14.181) . The scaling dimensions agree and the only difference is the 
eigenvector for 1/3. We have seen a similar phenomenon for the ^Sq channel, as we explained 
in Sec. HCl 
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By substituting 



\^6w J 



Ar 



A 



aui — + bu2 — 



Ao 



A 



A 



(4,41) 



into the amplitude, we obtain the off-shell amplitude near the nontrivial fixed point, 



27r2 V9 



2^ 



2b 



M , 



(4.42) 



where the ellipsis stands for higher order terms in a, b, and c. By comparing it with the 
effective range expansion. 



defined on the mass shell, one sees 



a 



M 
Air 



1 

a 



2A3 



:P + 



(4.43) 



8Ao^ 
" 97r 



(4.44) 



Namely, the two relevant directions correspond to the "scattering length" a and the "effective 
range" r. (These effective range parameters have different dimensionality in the P waves 
from those in the S waves.) Note that the strong coupling phase (a < 0) corresponds to a 
positive "scattering length," and the weak coupling phase (a > 0) corresponds to a negative 
one. 

In Ref. 30| it is noticed that there are two couplings to be fine-tuned. They corresponds 



to our relevant parameters given above. We have given the explanation for what they found 
from the RG point of view. 

In order to understand the difference between the strong and the wealc coupling phases, 
we take a closer look at the poles of the amplitude in the effective range expansion. Note 
that we are in the vicinity of the nontrivial fixed point, i.e^ \a\ -C 1 and \b\ <^ 1. At low 



energies, the poles are the solutions of the cubic equation 



30|, 



a 2 



0. 



Note that if we change the signs of a and r simultaneously, a 
solutions are given by p — > . 



-a and r 



(4.45) 
-r, the 
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Im[p] y 


c 








Re[p] 


/ 


\ 




X 


X 





Im[p] )^ 



X 



Re[p] 



FIG. 2: Poles for small < and r < 0. 
The arrows indicate the directions to which 
the poles move as 0—. There is a 

shallow resonance. 



FIG. 3: Poles for small > and r < 0. 
The arrows indicate the directions to which 
the poles move as 0+. There is a 

shallow bound state. 



Let us first suppose that we are in the region a ^ < but not very close to the critical 
value = 0, i.e., \ar^\ < 54. In this case, we have three solutions. 




1 



i I lap 
- r + 
6 \ 2v 



2|a|3 



1 



where 



V = 



108 + ar'' + 108a/ 1 



'54 



. (4.46) 



(4.47) 



is a positive number. These poles are in most cases unimportant in low-energy scattering. 
Only when we take the limit and r keeping the condition \ar^\ < 54, 

all of the three poles move closer to the origin and they become significant. The very 
existence of such a special limit is a consequence of the two relevant directions. Note 
that this limit is a simultaneous limit a ^ and b ^ with keeping \b^/a\ "small", i.e., 
\b^/a\ < 2 1 877r^/128 ~ 168.5. We do not consider such a case here, but approach to the 
critical value — 0— keeping r finite. This is the limit a ^ keeping b finite so that 
\b^/a\ oo. The discussion about these limits and the relations to the power counting 



suggested in Refs. (30|, l3l| is given in Sec. |Vl 
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If r < 0, the expression for the solutions of Eq. fl4.46p is still valid for \ar^\ > 54. The 
first two poles in Eq. (14.461) move toward the origin as we decrease the magnitude of 
and reach it at a^^ = (See Fig. [2]), while the third one does not move so much. This fairly 
deeply bound state is insignificant for low-energy scattering. Note that the first two poles 
represent a shallow resonance. 

If we go across the critical value = into the region where > 0, the two 
(degenerate) poles at the origin split on the imaginary axis. We have the solutions 



r 

6 



where 



2 sin 



Arg 



-r I 1 ± v3 cos , 
6 V V3 



sm 



(108 + ar^) i + 108 




(4.48) 



(4.49) 



The first pole represents a shallow bound state. See Fig. [3l The pole with the positive sign 
of the second corresponds to the bound state found in the < case. The pole with 
the negative sign of the second does not seem to have a definite physical meaning, but is 
responsible for the rapid change of the phase shift at low momenta. 

Note that there are several cases (the case with r > and a > and the case with r < 0, 
a > 0, and \ar^\ < 54) in which the poles appear on the upper half plane, but causality 



prohibits poles on the upper half plane except on the imaginary axis 42 



43 



44 



45|. One 



should there: 
Refs 



46 



47 



bre exclude such regions of a and b which give rise to acausal poles. See 



48 



49 



50 



5l[ | for the cases where some physical axioms such as causality, 



positivity and unitarity constrain the values of low-energy constants. Remember that we 
have assumed that \a\ <^ 1 and \b\ ^ 1 (and \b^/a\ is large), that is, <^ Aq and 

|r| <^ Aq, with \ar^\ being large, through the analysis to meet the condition that we are 
in the vicinity of the nontrivial fixed point. The poles we have discussed have magnitudes 
smaller than the cutoff Aq to be consistent with the effective range expansion. 

To summarize, we have done the RG analysis based on the cutoff independence of the 
amplitude, and found that, near the nontrivial fixed point, the scattering length and the 
effective range are related to the deviations a and b from the nontrivial fixed point. The 
amplitude has three poles at low momenta. Two of the three poles are sensitive to the value 
of and get close to the origin near the critical value. For r < 0, there is a shallow 
resonance in the weak coupling phase a > (a~^ < 0), and a shallow bound state in the 
strong coupling phase a < {a^^ > 0). The appearance of a (shallow) bound state in the 
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strong coupling phase is similar to that in the S waves. But the appearance of a (shallow) 
resonance is a new feature which is absent in the S waves. 



C. D waves 



There are four channels in the D waves; ^D2, ^-Di, ^D2, and ^D^-^Gs, the second of which 
mixes with ^Si and has been considered. We introduce the following projection operators, 



P. 



' 16^2 

(3^3) _\/T5 



ijk 



48 



V Vm Vi (icr2(T„)(zr2) 

1 



V i V J - I V ) (io-20-fe) + (cyclic in j, k) 
1 



g I y ^ ' ~ § ^"^^^^ ) + (cyclic in z, j, A;) [> (z(T2cri) 



PA 



7^5 



VjVjVfcV; — - V^( VjV j(5fc/ + (5 terms to be symmetric in i, j, k, I) 



+ ^ {Sij6ki + 5ik6ji + 6u5jk) 



{ia2(Ji){iT2). 



(4.50) 



Because the RG equations are the same for both ^D2 and ^D2 channels, we present the 
calculation for the ^D2 channel only. (Those for the coupled channel ^D-^-^G^ are similar to 
those for the ^Si-^Di and the ^P2-^F2 channels that we do not show them explicitly.) For 
this channel, the interaction terms for the Fa up to 0{p^) may be written as 



■^int 



+ 



+ 2B, 



(4,51) 



In terms of the dimensionless coupling constants. 



Xil,2) 



2n 



2 '-'4 



?/(l,2) 



27r2 ^6 



^(1,2) 



27r2 



r5 



(4.52) 
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we find the following RG equations (in the n — oo limit), 



^i^^ = - 5X(i,2) - [X(i_2) + 2X(i,2)y(l,2) + 2X(i,2)2;(l,2) + yfl,2) + '^V (1,2) Z (1,2) + 4,2)] 



dy(i,2) 
dt 

dZ(i,2) 

dt 



72/(1,2) - 
7^(1,2) + 



^a;(i,2) + 2a;(i,2) 2/(1,2) + ^2/(^2) + y{i,2)Z(i,2) - ^4,2) 
^a;(i,2) + X(i,2)y(i,2) - a;(i,2)2(i,2) + ^4,2) - y{i,2)Z(i,2) - ^4,2) 



(4.53) 



Just as for the ^Pi case, with the appropriate replacement of the coupling constants, the 
coefficients of the quadratic terms of the above RGEs are the same as those in the ^5*0 chan- 
nel, and the only difference is the coefficients of the linear terms representing the canonical 
dimensions. 

The nontrivial fixed points are found to be 



(^*l,2)5 2/*l,2)5 ^(1,2)) — (0;0;0), 



-5, 



25 25 



49 63 7 

'Y'To'~2 



(4.54) 



6 6 

but the last nontrivial fixed point has complex scaling dimensions, and thus may be disre- 



garded. At the nontrivial fixed point, (— 5, ^, — ^) , we find the following scaling dimensions 
and the corresponding eigenvectors, 

/ 3 \ / n \ 



/-lo^ 





ui = 5 : Ui = -5 , z/2 = 3 : M2 = -1 , = -2 : = 5 . (4.55) 

V W V 2 / 

Note that there are two relevant couphngs. For the first two eigenvectors, the scaling di- 
mensions shift from their canonical values by ten. 

To summarize, we have examined the D waves in the similar way, and found that there 
are two relevant couplings with the scaling dimensions shifted by ten. These large anomalous 
dimensions suggest that there may be more relevant operators. In the next section, we give 
another guide for what we expect for scaling dimensions from a different point of view, and 
show that there is one more relevant operator. 



D. PDS for higher partial waves 



In the previous paper [18|, we explain how the PDS leads to the "shift by two" rule for 
the S waves, that is, the anomalous dimensions of the four nucleon operators are two. In 
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short, the usual PDS renormahzation, which subtracts the contribution at the D = 3 pole 
as well, treats the operators as if they were in three dimensional spacetime. The canonical 
dimensions of the four-nucleon operators shift by two, e.g., the operator {N'^Ny{N'^N) has 
dimension six in (1 + 3) dimensions, but four in (1 + 2) dimensions. (In D dimensions, it 
has dimension 2{D — 1).) 

Note that this "shift rule" does not apply to the redundant operators, which do not need 
to be introduced with the dimensional regularization. 

We have seen that the rule for the P waves is "shift by six," and for the D waves, "shift 
by ten." It is now natural to extend the PDS scheme to higher partial waves by subtracting 
the contribution at D = 1 for the P waves, and that aX D = —1 for the D waves. 

One can easily show that this generalization of PDS for higher partial waves works well 
as for the S waves in the pionless NEFT. 

This kind of generalization of PDS has been considered in Ref. 5^ but apparently the 
relevance to higher partial waves was not noticed. 

Guided by this extended PDS prescription, we expect that there should be one more 
relevant operator (with the scaling dimension one) in the D waves, even though we found 
only two relevant operators in the previous subsection. The operator will be found if one 
performs a similar analysis in larger operator space. 

Note that the RGEs for the P and D waves have the similar structure to those for the S 
waves, as we remarked in Sec. IIV Al and in Sec. IIV Ci The RGEs for the D waves to 0{p^) 
should be obtained from those for the S waves to 0{p'^) given in Sec. IIIIBI by replacing 
the coefficients of the linear terms with the corresponding ones for the D waves (which 
are nothing but the canonical dimensions), with appropriate replacement of the coupling 
constants. By examining the RGEs, one can immediately obtain the nontrivial fixed point, 

25 25 100 100 100\ 
-5,-,--,0,— ,0,- — ,-— J, (4.56) 

where the first three components are (a:(i,2), y{i,2), ^(1,2)), while the rest are the dimensionless 
coupling constants of 0{p^). 

By solving the linearized RGEs around it, we obtain the following eigenvalues and eigen- 
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vectors, 

J^i = +5: ni = 0,-1,1,0,-1,0,1,1^ , (4.57) 

zy^ = +3 : U2= (^0, 0, -1, 0, 1, 1^ , (4.58) 

z/3 = +1 : = (0, 0, 0, 0, -1, 0, 1, 1) , (4.59) 

/ 54 90 261 190 64 \ 

''' = -^'- =1-127' 127' -635'°' -127'°' 127' V' ^ '^^ 

/ 30 25 25 27 \ , , 

-. = -4: -2. 0.0, ij. (4.61) 

..^-4: 1.0.1.0). (4.62) 

/ 45 55 55 155 \ ^ 

^.--4: ...^(-^,-^,^,^,0,l,0,0j, (4.63) 

= -9 ; tig = f -, --, A, 0, -, 0, -, l] . (4.64) 

^ ^ V8 8 10' ' 8 ' ' 50' y ^ ' 

We find the third relevant operator as we expected. The extended PDS prescription and the 
Wilsonian RG analysis given above give a consistent result. 



V. SUMMARY AND DISCUSSION 



In this paper, we extend our previous study of the Wilsonian RG analysis for the pion- 
less NEFT in the two-nucleon sector. The determination of power counting based on the 
scaling dimensions is a simple and powerful method in particular for the theories in which 
nonperturbative dynamics is important. 

Two kinds of extensions are considered; (1) enlargement of the space of operators to be 
taken into account, and (2) higher partial waves. Because our formulation is general, we can 
use the same machinery to analyze them. 

We considered the space of operators up to including those of C(p^) in the ^Sq channel 
and found that the results are stable against the enlargement. We also found that there are 
"fixed lines" and a "fixed surface" which are related to the existence of marginal operators. 

In the P and D waves, we derived the RG equations and found the phase structures. 
There are two phases, the strong coupling and the weak coupling phases, just as in the 5" 
waves. Unlike the case of the S waves, however, there are two relevant directions at the 
nontrivial fixed point for the P waves. We found that there are three for the D waves. 
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By explicitly calculating the off-shell amplitude for the P waves, we have seen that (near 
the critical surface and r < 0) there is a shallow bound state in the strong coupling phase, 
while in the weak coupling phase there is a shallow resonance. 

To summarize, we have a coherent picture of the pionless NEFT in the two-nucleon sector 
from the Wilsonian RG point of view. 

In the following, we discuss several aspects of the results. 

1. In the enlarged space calculation, we found "fixed lines" and a "fixed surface." At 
first sight, they seem very strange. But actually, their existence is related to that of 
marginal operators. In relativistic field theory, classically marginal operators usually 
get (non-integer) anomalous dimensions and turn into relevant or irrelevant operators, 
so that the corresponding coupling constants run though slowly. In the present case, 
we have no logarithmic divergences so that the marginal operators are really marginal. 
Their couplings do not run at all. 

2. The absence of logarithmic divergences seems to lead to integer scaling dimensions. 
We found however that "fixed surface" [F] has irrational scaling dimensions. This fact 
tempts us to think that it is an artefact. 



3. The additional nontrivial fixed point found in Ref. [18|, which has complex scaling 
dimensions, disappears in the enlarged space calculation. It is well known that the 



3- 



It is also known that trun- 



truncation in general produces spurious solutions [3i 
cation leads to complex scaling dimensions. (Note that the RG equations obtained 
from amplitudes do not have such spurious fixed points.) As we discussed in Sec. IIII Bt 
those fixed points with complex scaling dimensions are considered as spurious. Note 
however that what we are doing is not the so-called polynomial expansion (the expan- 
sion in powers of fields), but the derivative expansion only. It is due to the Fermi-Dirac 
statistics and the nonrelativistic feature that there are only a finite number of inter- 
action terms in the derivative expansion to a given order. It is remarkable that the 
derivative expansion shows the similar symptom to that of the polynomial expansion 
seen in relativistic scalar theory. 



4. As we discussed in a previous paper [18|], some of the eigenvectors seem to correspond 
to the directions in which the physical quantities remain unchanged. We suspect that 
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the eigenvectors with scahng dimensions which do not obey the "shift by two" rule 
in the S waves are such directions, and that there is a similar correspondence in each 
partial wave. 

5. In the P wave calculation, we found that there are two relevant operators in a very 
natural way. The fact that there are two couplings to be fine-tuned has been noticed 



in Ref. [30|]. It requires a nonperturbative analysis to decide which couplings should 
be treated nonperturbatively and cannot be determined by perturbative considera- 
tion. In fact, they showed that two couplings should be fine-tuned by the explicit 
(nonperturbative) calculation using a dimeron. 



6. In Ref. [30t|, the case 

^ Ml, r ^ Mio, so that lar^l ~ C(l) (5.1) 



is considered as an unnatural case, where Mio is a low-energy scale. In Ref. 3l|, on 
the other hand, another power counting is considered, in which only one combination 
of coupling constants is fine-tuned. 



a 



-^r^MlMh^, t^Mm, so that jar^l ~ O (^^J > 1 (5.2) 



where M^i is a high-energy scale. The former case corresponds to the special limit 
mentioned in Sec. lIVBt in which is sent to zero, keeping |ar^| < 54. All of 
the three poles move closer to the origin, thus they become significant in low-energy 
scattering. The latter corresponds to the limit in which — > keeping r finite, 
so that |ar^| > 54. Only two poles (representing a resonance in the weak coupling 
phase) move closer to the origin, while the third is insignificant. In terms of the phase 
diagram, this corresponds to the flows very close to the phase boundary, but not to 
the fixed point itself. 

7. In the real world, the two-nucleon system in the P waves does not exhibit any shallow 
resonances nor bound states, so that it is in the weak coupling phase. There may 

i^T at low energies. For example. 



be systems which can be described by the same E 

Feshbach resonances of ultracold {5^ or ^Li 5^, |55| may be interesting 
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APPENDIX A: RELATION BETWEEN THE LEGENDRE FLOW EQUATION 
IN THE SHARP CUTOFF LIMIT AND THE RG EQUATION BY BIRSE ET AL. 



1. Feynman rules 

Legendre flow equation reduces to a set of RG equations for coupling constants, which 
consists of one-loop diagrams. In the present case, a typical diagram is given in Fig. |H Let 
us first describe the Feynman rules for the ^5*0 channel as an example. From the averaged 
action, one can easily read off the Feynman rules for the vertex (See Fig. [5]), 

AtXAFAiPf,P^){Pl),,iPa\,, (Al) 

where xa is a dimensionless coupling constant. For the ^5*0 channel, it is one of the dimen- 
sionless coupling constants introduced in (13. 7p . xa G {x,y, z,ui,U2, zi, Z2, z^}. {Pa)ij is the 
spin-isospin factor of the projection operator to the partial wave in question, with the indices 
i and j referring to the spin and isospin quantum numbers of the nucleon pair. FA{pf,Pi) 
is the corresponding momentum-dependent factor to Xa, where Pi stands for the incoming 
momenta, while pf for the outgoing ones. For the ^5*0 channel, we have 

i=l 

-27r^ / 2 2 \ p _ -27r^ 
^""^ ~ 16MA^ ^""^^ + ""^^^ ' " TgMA^'"''''''' 

Pzi = ^5 ^ Si, = ^5 (5*1 + 5*2) {S3 + 5*4) , 

i=l 

F., = ^ {ri2 (^3 + ^4) + r34 (^1 + ^2)} . (A2) 
Note that r^j and Si are defined in Eq. (13. 6p . 



2. Subtleties at equal times 

When we include higher order redundant operators, we encounter subtleties which do not 
emerge in the lowest order. In this subsection, we explain the subtleties and how to handle 
them. 
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I Pi 



j P2 



FIG. 4: One-loop diagram contributing to 
the Legendre flow equation. Labels A and 
B stand for the type of the vertices. The A 
vertex is given in FIG. [5j 



FIG. 5: Feynman rule for a vertex labeled 
by A, where i, j, etc. are the spin-isospin 
indices, and {k = 1,2,3,4) are momenta 
of the nucleons. 



The Legendre flow equation may be obtained by considering the following integral 
\{Aixa) (4zxb) [Pi) ^^Ti [PaPl] (Pi 



b)ij 



pO/2 + A;0-7^((p,/2 + fc)') 



le 



le 

Fsik^pf), (A3) 



where we have introduced 



and 1/2 is the symmetric factor. (Note that the way the IR cutoff functions appear should 
be taken as symbolic, because it comes from the "naive" inclusion of the cutoff function in 
the averaged action that however breaks Galilean invariance. A precise, Galilean invariant 
way is given shortly.) 

The problem is that the product of the factors Fa and Fb may be quadratic or higher 
order in k^, so that the integral over appears to diverge. This divergence cannot be 
regularized by the introduction of the higher order terms in k^ to the denominator of the 
propagator, because the results depend on the inverses of (expectedly small) coefficients. It 
means that the higher order terms drastically change the lower order results, and cannot be 
accepted. 

Note that in the relativistic field theory with dimensional regularization such divergence 
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does not cause a problem, because of analytic continuation. 

The divergence comes from the large region, in other words, from the infinitely short 
time intervals. But from the EFT point of view, one should have a cutoff on the energies 
of the intermediate states. (There should be a finite resolution of time.) The cutoff may 
be of order Aq/M, where Aq is the physical cutoff on three-momenta. We therefore assume, 
though implicitly, that is cutoff at a scale of order Aq/M and ignore the divergence arising 
from the k^ integration. Effectively, it ends up with evaluating the integrand at the (either) 
pole. 

The manipulation of making the cutoff Galilean invariant is evident at the amplitude level: 
to impose the cutoff on the relative momentum [l^. After doing so, Eq. (1A3P becomes 

where we have introduced 

and FA{pi, k)FB{k,pf)\p^i^ denotes that Fa^Pi, k)FB{k,pf) is evaluated at the pole. 
3. Equivalence 

In order to establish the relation between the Legendre flow equation and the RG equation 
employed by Birse et al., it is useful to rewrite the Legendre flow equation in the sharp cutoff 
limit in a simpler form. 

In the sharp cutoff limit, the integral in Eq. ( 1A5|) may be further simplified, and can be 
written as 

A- (Pi\ (P\ 1 /•.20 FA{p^,k)FB{k,Pf)\^^i^ 

- 4^XAXs (Pi),, {Pah, —,J^ kdkjd , (AT) 

where we have used Tr (PaPb) = ^abf^- The change of A may be compensated by the change 
of the coupling constants. Thus, if Cc ~ xc/A.'^^ , 



dt 

receives the contribution from 



^^'^ Fcipf,Pi) + dcXcFcipf,Pi) (AJ 



in-" 



^ \ dAj J " A{pi) - k^/M + le 



MA f ^a(p.,A)Fb(A,P/)|,,,, 
Stt^ i 1 - A{pi) 
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where A{P) = MA{P)/h? and ^^(pt, A)Fb(A, p/)!^^^^ is the simphfied expression for 
^AiPi, k)FB{k,pf)\^^^^ with the magnitude of k being A. In the center-of-mass frame, all 
the Fa in the ^5*0 channel do not depend on angles, so that we obtain 



dxc , , MAFAip^,A)FB{A,Pf\,l, 



dt ^ 27r2 l~A{j)i 



(AlO) 



c 

where \c stands for the operation of taking coefficient of Fc{pf,Pi) in the expansion of the 
right hand side. This is a very handy expression for the sharp cutoff limit. 
If we identify the "potential" V{pf,pi] A) as 

Vipf,p,;A) = -J2^cFc{Pf,Pi), (All) 
c 

we have 



dV _ 
aA ~ A 



c ^ ^ c 

Substituting Eq. ( lAlOl) . we have 



dV M ^ ^ . ^, 1 ^ 



A ^ B 



=^^np,,A;A) ^,J,^^^/ (A,KAl (A13) 

Since A{pi) = p° in the center-of-mass frame, it is nothing but the RG equation (12. 3p 
employed by Birse et al. 



APPENDIX B: CUTOFF FUNCTION DEPENDENCE OF THE RESULTS 

In this section, we briefly show how the results for the ^5*0 channel depend on the param- 
eter n in the cutoff function R\ in Eq. (12. 2p . We have derived RG equations for an arbitrary 
value of n, but the expressions look too complicated that we omit them. In the following, 
we refer the fixed points/lines/surface as [A], [B], etc., as we defined in Sec. IIII B[ 

The primary objective of this study is to justify the use of the sharp cutoff limit, which 
gives the simplest expressions for the RG equations. It is known that the sharp cutoff leads 



to bad behaviors in the derivative expansion 36[], so that it is important to see how the 



results behave as we approach to the sharp cutoff limit. In particular, we are concerned 



34 



TABLE I: The dependence on n of the fixed point [B]. 



n location of the fixed point scaling dimensions 



2 




1 in'^9fi 

J- . lUO^U , 




— 1 .U I U, -L .O / UU^, 1 .U I UU^i J 


(' — 5 


—4 


—4 


—4 


•J) 


_9 


-1 1) 


10 


{- 


1.02722, - 


-0.524991,0.524991,0, 


-1.40602, 0, 1.40602, 1.40602) 


(-5, 


-4, 


-4, 


-4, 


-3, 


-2, 


-1,1) 


102 


(- 


1.00287, - 


-0.502587,0.502587,0, 


-1.34072, 0, 1.34072, 1.34072) 


(-5, 


-4, 


-4, 


-4, 


-3, 


-2, 


-1,1) 


103 


( 


-1.00029. 


-0.50026,0.50026,0,- 


1.33407, 0, 1.33407, 1.33407) 


(-5, 


-4, 


-4, 


-4, 


-3, 


-2, 


-1,1) 


10^ 


(- 


1.00003, - 


-0.500026,0.500026,0, 


-1.33341, 0, 1.33341, 1.33341) 


(-5, 


-4, 


-4, 


-4, 


-3, 


-2, 


-1,1) 


CO 






\ 2 ' 2 ' ' 


3j ^> 3) 3; 


(-5, 


-4, 


-4, 


-4, 


-3, 


-2, 


-1,1) 



with the possibility that the sharp cutoff may cause non-analytic, or singular behaviors. 
We therefore consider the smooth (i.e., differentiable) cutoffs which are very close to the 
sharp cutoff limit. At any rate, we do not pretend to settle the convergence problem caused 
by sharp cutoffs in the expansion. We only demonstrate numerically that the use of sharp 
cutoff does not seem to cause a serious problem in our present case. On the other hand, some 
authors try to extract the information about the convergence of the derivative expansion by 
looking at the so-called "scheme (in-)dependence" 29|]. From our point of view, the "scheme 
independence" is a necessary, but not sufficient condition for the convergence. The best we 
can do is to actually enlarge the space of operators as we are doing in this paper. 

The trivial fixed point [A] is always there, and has no n dependence at all. The n 
dependence of the nontrivial fixed points [B] and [C] is given in Tables [T] and [iTl The fixed 
point [B] is very stable against the variation of n. Even though the location changes slightly, 
the scaling dimensions do not change at all. The fixed point [C] has a stable limit but it has 
complex scaling dimensions. 

We are not completely sure that the "fixed lines" and the "fixed surface" exist for an 
arbitrary n, because we are unable to obtain the analytic expressions for them. Numerical 
study indicates that at least "fixed line" [E] is stable. In Table IIIIl we present the n 
dependence of a particular point {u2 = 0) on the line and the scaling dimensions at the 
point. 
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TABLE II: The dependence on n of the fixed point [C]. means that scahng dimension i/j is 
doubly degenerate. 



n location of the fixed point scahng dimensions 

2 (-2.03159, 1.46753, 1.27909, 0., -2.64366, 0., 1.78644, -1.84813) (-6.36576, (-5.68288)2, 5, 3, 1, -1.68288 ± 2.27988i) 

10 (-2.27317, 2.69046, 1.2095, 0., -5.60121, 0., 1.3641, -2.62099) (-7.85177, (-6.42588)^, 5, 3, 1, -2.42588 ± 3.12036i) 

10^ (-2.256, 2.79556, 0.984983, 0., -5.98406, 0., 0.934766, -2.2433) (-7.99819, (-6.49909)^, 5, 3, 1, -2.49909 ± 3.20057i) 



10'^ (-2.25064, 2.79559, 0.957635, 0., -5.99646, 0., 0.888947, -2.19038) (-7.99998, (-6.49999)^, 5, 3, 1, -2.49999 ± 3.20155i) 
10"^ (-2.25006, 2.79547, 0.954855, 0., -5.99738, 0., 0.884356, -2.18494) (-8, (-6.5)2, 5^ 3^ _2.5 ± 3.20156i) 

i-l 0^ -fii^ 0^ iH. -Hi) (-8, (-13/2)^ 5, 3, 1, (-5 ± ^V4^)/2) 

TABLE III: The dependence on n of the fixed point ii2 = on the line [E]. uf means that scaling 
dimension z/j is triply degenerate. 



n location of the fixed point scahng dimensions 



2 


(0,0,0,0,0,0,4.41305,0) 


((-5)^ 


-3, 


-1,0,2,5) 


10 


(0,0,0,0,0,0,5.51631,0) 


((-5)^ 


-3, 


-1,0,2,5) 


102 


(0,0,0,0,0,0,5.07010,0) 


((-5)^ 


-3, 


-1,0,2,5) 


10^ 


(0,0,0,0,0,0,5.00719,0) 


((-5)^ 


-3, 


-1,0,2,5) 


10^ 


(0,0,0,0,0,0,5.00072,0) 


((-5)^ 


-3, 


-1,0,2,5) 


CX3 


(0,0,0,0,0,0,5,0) 


((-5)^ 


-3, 


-1,0,2,5) 
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